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ABSTRACT 


This thesis describes a circuit model for the inductive strip 
in inhomogeneous finline with the following geometry: fin and strip 


centered in the shield, dielectric material with €,=2.22, 


b/a =4/9, 0.5 S W/b S$ 1.0, T/a 2 0.01 and 0.0 $ d/a S$ 0.1 . The 


model is shown to produce results that agree with data computed 


using the spectral domain method. The model has been generated 


using WR(90) waveguide operating in the TE,, mode. 
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I. INTRODUCTION 


A. BACKGROUND 

Finline is a transmission structure for electromagnetic 
waves that was first discussed by Meier in 1974 [Ref. 1]. 
Finline consists of one or more thin metal fins printed on a 
dielectric substrate mounted in the E-plane of a rectangular 
waveguide. Figure 1 depicts the particular variation under 
consideration here. There are several advantages to the 
finline structure. The fin manufacture is simplified by the 
presence of the dielectric substrate which allows the use of 
well developed etching technologies for their manufacture. 


Finline has less_ stringent tolerance requirements’ than 


Inductive Strip in Finline 


= Metal Fin and Strip 





Figure 1. Arrangement of inductive strip in finline. 








Microstrip and the use of a rectangular waveguide shield 
simplifies the interfaces with other rectangular waveguides. 
Like conventional waveguides, finline can operate in single 
mode with low attenuation. [Ref. 1] 

Computer assisted design (CAD) is essential to the 
development of complex devices in all fields of engineering. 
This is especially true in microwave and millimeter wave 
design. Numerical solutions to many electromagnetic problems 
result in time consuming computer implementations. To be 
useful in design, a workstation should be able to quickly 
predict the response of a known structure. The use of 
equivalent circuit models for common structures that respond 
in the same way as the general solution for a limited, but 
useful, range of parameters has been found to be a useful CAD 
technique. 


B. USE OF PREVIOUS WORK 

Initial work by Knorr and Shayda [Ref. 2] solved the 
electromagnetic fields in an arbitrary section of finline 
using the spectral domain method. The program IMPED 
implemented the solution. Knorr and Deal solved for the 
scattering parameters of an arbitrarily located inductive 
strip in finline and implemented the solution in a program 
called STRIP [Ref. 3], (Ref. 4]. Morua developed a circuit 
model for homogeneous finline, using input from both IMPED and 
STRIP {Ref. 5). Grohsmeyer developed a model for 
inhomogeneous finline using the IMPED program. This work uses 
all of the previous work as a basis to begin and relies on the 
work of Grohsmeyer for the finline model [Ref. 6]. 


C. OBJECTIVE 

The purpose of this study is to create a model that allows 
the representation of an inductive strip in inhomogeneous 
finline. The model must be accurate and suitable for 


implementation in a CAD environment. 








The desired range of validity is as follows: 


°* €& = 2.22 

e W/b - ratios from 0.5 to 1.0 

°e b/a - ratio of 4/9 

¢ T/a - ratios greater than 0.01 
e d/a - ratios from 0.0 to 0.1 


D. CONVENTIONS 

Previous work, referenced above, specified the subscript 
‘eq’ to denote an equivalent dimension (length or dielectric 
constant). This model requires several equivalent dimensions 
to be specified. Therefore, the following subscript 


convention will be used. 


e ’ ' = actual dimension (no subscript) 

e ‘f’' — equivalent dimension in finline 

e ‘dad’ — equivalent dimension in dielectric loaded half 
waveguide 

e ‘a’ - dimension of the air filled half waveguide. 


Figure 2 depicts three views of the strip in finline and 
labels the parameters used to describe the geometry. The 
dotted lines are grid lines to assist in specifying the 
dimensions. 





Finline Dimensions 


Top View 


End View 


Side View 





Figure 2. View of the inductive strip in finline with 
parameter dimensions indicated. 
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height of the shield 
length of inductive strip 


° 
axqn Oo w 
| 


thickness of dielectric 


II. MODEL DEVELOPMENT 


A. CONCEPT 

The model for an inductive strip was derived from the 
homogeneous inductive strip model develoned by Morua [Ref. 5]. 
The model retains a similar structure but the calculation of 
element values has been modified to account for the inclusion 
of dielectric. 

The model takes the physical structure of the finline 
inductive strip and divides it into two half-waveguides that 


are treated as separate elements, as shown in Figure 3. One 


Two Waveguides 





Figure 3. Looking into each of the half waveguides. 





waveguide is air filled with transverse dimension a/2 and b. 
The second waveguide is dielectric loaded with the same 
dimensions as the air filled guide. At the mouth of each of 
these half-waveguides, there is a parallel inductor/capacitor 
circuit that accounts for stored electric and magnetic energy. 
The circuits are then connected through a perfect impedance 
transformer to the finline. Figure 4 shows the entire 


structure. 





Figure 4. This figure shows the circuit at one end of the two 
half-waveguides. The other ends of the waveguides are 
connected to the mirror image of the above circuit. 


B. DERIVATION OF THE MODEL 
1. General Comments 
A model for the loaded half-waveguide that could be 
represented in CAD software was found. This required that an 
equivalent below-cutoff homogeneous waveguide be found. A 
model form that replicated the scattering properties of the 
inductive strip as predicted by STRIP was then developed. 
2. Loaded Waveguide Model 
The loaded, below-cutoff waveguide was modeled by 


matching the propagation constant and voltage-power impedance 





with that of equivalent homogeneous waveguide. The first step 
of this process was to write the solution for the 
electromagnetic wave propagation for a below-cutoff dielectric 
leaded waveguide as shown in Figure 3. In the following 
discussion, it is understood that, B, will be a positive 
imaginary number due to the below-cutoff operation of the 
guide and will result in a negative real propagation constant. 
The full derivation of this solution is shown in Appendix A 
and parallels the discussion on propagating loaded waveguide 
in Pozar, [Ref 7;pp. 151-153]. This analysis results in the 
following transcendental equation 


o=/k¢-B? tan (dfe,ke-B? ) +/e,ke-B? tan((a-d)/k7-B? ). (1) 


Since this equation must be true for all frequencies, it is 
possible to find a set of, $’s, for various frequencies and 
geometries. Several MATLAB functions were written to compute 
B, in a half section of WR(90) waveguide for frequencies 
between 6 GHz and the cutoff frequency, f,, and for 2d/a 
values from 0 to 0.25 . Figure 5 shows the cutoff frequency 
as a function of d/a for WR(90) waveguide. Figure 6 depicts 
B versus frequency for various dielectric thicknesses. These 
curves describe the propagation characteristics of a 
particular geometry and are the characteristics that need to 
be matched by the equivalent waveguide. The equation for B in 


a homogeneous waveguide can be written as 


B= lens {2 (2) 
Gq 


Using another set of MATLAB routines values of a, and €,, that 
cause a waveguide described by equation (2) to match the 


curves in Figure 6 can be found for each desired geometry. 
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Frequency 


GHz 


in 





frequency for several values 


vs. 





That data can then be reduced to the following two equations 


using a least mean square error curve fitting routine: 


€ =1.0+0.0192 (2-2) -0.3391 (2-2) 2410.6174 (2-2) 34 
rd a a a (3) 
-3 6647 (2-2) 4+0.01356 (22)5 


and 

a 

24 21 0+0,001314 (2-2) -0.03026 (22) 240.3424 (2-2) 3+ 

a/2 a a a (4) 

-2.1617 (29) 441 ,.7461(22)5 
a a 
where 
oxs2<0.1. 
a 


The shapes of the curves described by the above equations are 
shown in Figures 7 and 8. Application of these equivalent 
dimensions to equation (2) results in B’s that differ from 
those found for the loaded waveguide by less than 0.2%. 
MATLAB functions used in the production of the equivalent 
homogenous half-waveguide in this section can be found 
Appendix B. 

The voltage-power impedance of the loaded waveguide 
can be calculated analytically since the electric and magnetic 
fields can be found to within a single constant value. The 


voltage-power impedance is given by 


Zev. (5) 








€. Variation with 2d/a 
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Figure 7. Variation in €, with 2d/a. 
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Variation of a with 2d/a. 


Figure 8. 
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In this case, since the waveguides are operating beiow-cutoff 
the power will be imaginary and will result in an imaginary 


impedance [Ref. 8]. The power can be found to be 


pa Aeros tan(k,(a-d) | (a-d) A?wBp 
2k; 2kzcos (k,(a-d) ) 
: ° ; (6) 
_ A? woBpsin (kd) cos (k,d) , A?dw By 
2k3 2k3 | 


Equation (6) can be written more compactly as 
P=b C, (7) 


where the imaginary constant C, represents the terms in 


parentheses above. The voltage in equation (5) is given by 


b 
v= -[E, d1=-E,b (8) 
{?) 


where E, is the maximum electric field. The voltage-power 
impedance for the loaded half-wave guide is 
2 


E 
Zoy=D— : (9) 


The voltage-power impedance for the equivalent below-cutoff 
waveguide is given by [Ref. 6] 


2bg ‘ 
ag/E-g ae x 2 (10) 
(zeae) 





Zoy2120R 
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b, can be found by 


Given that a, and €,, are already known, 


for b; and substituting in 2Z,, 


from 


(10) 
equation (9) which results in 


solving equation 





(11) 


and (11) define the loaded half-waveguide 


Equations (3), (4) 


The variation in by, with dielectric thickness is shown 


model. 


in Figure 9. 
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3. Integrated Model 

Using the circuit model developed by Morua as a 
starting point, a MATLAB function that replicated the 
performance of the CAD workstation was written. This program 
was validated using arbitrary values of the inductance, 
capacitance, and turns ratio on the CAD workstation. The CAD 
program and the MATLAB functions used are located in Appendix 
Cc. A function minimization routine was then written which 
found values of the various elements that resulted in the 
model output matching the spectral domain data produced by 
STRIP. During the course of much trial and error, it was 
found that the expression used by Morua to calculate the turns 
ratio in the homogeneous case could not be extended to the 
inhomogeneous case. The model produced good results when the 
conductance on the strip side of the transformer was specified 
as a function of frequency. It was found that the conductance 
could be defined in the following way 





Y)=0.001(t,f,+¢,f4+t,f£,+C) (12) 
where 
£=£ 
£5 (13) 


Figure 10 depicts the way conductance varies with frequency 
for a typical geometry. Once the value of the conductance 
needed to match the spectral domain data was found, the turns 


ratio of the impedance transformer was calculated using the 


relation 











where N is the turns ratio, Y, is the finline conductance, Y, 
is the conductance on the strip side of the transformer. The 
behavior of the turns ratio with frequency is shown in Figure 
11. This model was found to be able to replicate the 
performance of the inductive strip for all geometries of 
interest with less than 1% error. Some results are shown in 
Figures 12 and 13 (in the figures spectral domain data is 
referred to as SPEC_DOM and the model is referred to as 
STRIP). 


Conductance Variation with Frequency 
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Figure 10. Variation of the conductance with frequency. 











ry 
¢c 
© 
bm] 
ica 
2 
ke 
<= 
= 
= 
2 
~~ 
9 
m 
a 
© 
e 
-] 
i 
£ 
= 
2 
~~ 
2] 
8 
> 


ee en Se Se eee 


mee mee et ete tee eee pen ee ee ee qe nnn ene ene genre cnn sce epee n ence eens ene 


eh ee bh ne nn te Nn ne ee bee ee ee bee eee ee new eee 


” 
Nn 


ooy Suny 


Nee ee eee eee ere eee ee eee eee ee eee eee 


W/b=0.5, T=50, and d=350 


Frequency (GHz) 





Figure 11. Variation in Turns Ratio with Frequency. 
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Figure 12. Response of S,, for Model and Spectral domain data, 
for W/b=0.5, T#=10 mils, d=0 mils, L=14.67nH, C=0.0043pF and 
Y,=0.001(1.315+3.63f,-3. 43f,7+0.98f,”) . 
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and Y,#0.001(1.69+2.06f,-1.52f,,+0.3792f,°) . 
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III. DATA COLLECTION 


A. STRIP PROGRAM DESCRIPTION 

The input data was generated using the STRIP program. 
STRIP is a Fortran program that computes the scattering 
parameters of the inductive strip in finline using the 
spectral domain method. The STRIP program and the development 
of the theory behind it are discussed in detail in References 
3 and 4. Accuracy of 1-2% with experimental data on 
homogeneous finline were reported. [Ref. 3], (Ref. 4] 


B. MODIFICATIONS TO THE ORIGINAL PROGRAM 

During the course of the data collection both a SUN 
workstation and a VAX 11/785 were available to compute the 
scattering data. The two computers were found to produce 
results that differed by as much as 10%. These errors were 
presumed to result from differences in the way the two 
machines handled small numbers. This indicated a possible 
problem in the way the program dealt with small numbers. To 
correct the problem, the program was changed to re-specify 
some variables that were previously real as double precision. 
This modification elevated the problem and resulted in 
agreement between the solutions generated by the VAX and the 
SUN. Unfortunately, other intermittent problems began to 
appear in the output data. Further extensive trouble shooting 
revealed that some double precision values were being returned 
from subroutines into real numbers. The extra bits were 
written into the adjacent memory spaces, corrupting them. 
This problem was corrected by re-compiling the program to 
treat all real numbers as double precision. This final 
alteration corrected all the apparent problems in the 


operation of the program. 
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C. DATA GATHERING 


The data on which the model was based, was calculated 


using parameters associated with WR(90) waveguide. WR(90) 
waveguide was selected for two reasons. First, it is more 
convenient to conduct measurements in X-band. Second, as 


shown by Morua in the homogeneous case, selection of b/a=4/9 
May permit the model to be extended to b/a’s, from 0.4 to 0.5. 
The following list specifies the values of the parameters used 
in determining the data points. All measurements in mils 

(1 mil = 0.001 inches): 


e W/b = { 0.5, 0.7, 0.9, 1.0 } 


°* T= { 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 
400, 500, 600, 700, 800, 1000 } 


edz= { 0.0, 10, 30, 50, 70, 90} 


The strip length interval was increased for T>100 mils because 
the scattering parameters were found to vary more slowly at 
large values of T. 

The performance of the STRIP program is determined in 
large part by the order of the matrix of inner product terms. 
Morua showed that most of the convergence for typical 
parameter values had occurred by matrix of order 10. 
Therefore, a matrix of order 10 was used as a starting point 
to calculate the data. Some data sets were re-computed with 
matrix orders of 12 to attempt to correct errors discussed 
below. 

The STRIP program was developed based on the assumption 
that the finline was operated in the TE, mode and the two 
half-waveguide sections were operated below cutoff. As the 
thickness of the dielectric substrate was increased the second 
propagating mode began to have an effect for the higher 
frequencies in the range of interest. In this situation, the 
program returns an increasing S,, magnitude normally 1.0. 


Values of the magnitude of S,, were assumed to be continually 
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decreasing with increasing frequency. This assumption 
complies with intuition of the physical system that as the 
non-propagating half-waveguide approaches the conditions for 
propagation, more energy will be transmitted and less 
reflected. Therefore, data points that indicated an 
increasing magnitude of S,, with frequency were discarded. At 
apparently random geometries, the STRIP program returned 
values for the magnitude of S,, that did not correspond with 
the similar, but different, geometries or the expected 
behavior of the structure. Identification of this data is 
discussed below. The block of erroneous data would typically 
begin at a strip length, T, of approximately 500 mils and 
continue through 1000 mils for a given dielectric thickness 
and gap width. This type of error was found to effect 15% of 
the data. By not including data with T greater than 500 mils 
the amount of affected data was reduced to approximately 5%. 
For this reason, data with a strip length greater than 

500 mils was discarded. An example of this behavior is shown 
in Table 1 (on the following page). The bad data is located 
at d=70 mils, T=600 mils and f=12 GHz. The numbers in the 
table represent the magnitude of S,,. The physical arrangement 
of the system is such that with increasing strip length, T, 
the reflection coefficient should continually increase. The 
adjacent data, at d=50, does exhibit the expected behavior. 
The data point described here and others that could be 


identified as bad in this manner were discarded. 











Table 1. Computed Magnitude of S,, near d=70 and T=600 
at 12 GHz, Illustrating Anomalous Results 


£=12GHz a=50 mils d=70 mils 
T=300 mils 0.8790 0.7877 
T=400 mils 0.9643 0.8216 





T=500 mils 0.9788 0.8242 
re6oo mits | 0.9872 | 0.6259 








III. ELEMENT REDUCTION 


A. CONCEPT 


It has been shown that element values could be found which 
produced excellent agreement with the STRIP data for specific 
geometries. The development of a useful tool for the CAD 
environment requires that a method of generating element 
values as a function of the geometry be found. In the 
following discussions, the coefficients of the strip 


conductance are referred to as elements. 


B. REDUCTION TECHNIQUE 


Polynomials were selected as generating equations. 
Polynomials were selected because of their simplicity and 
because the element data appeared to be slowly varying. The 
function minimization routines resulted in solutions that were 
exact beyond the level needed for a successful model. Errors 
of less than 1% with the STRIP data were typical and the 
expected error for the STRIP data was 1~2%. Therefore, it was 
assumed that an acceptable solution could be found and that a 
least mean square curve fit would be an adequate method of 
reducing the data. Programs were written that conducted a 
least mean square fit over each of the three parameters of 
interest. 

The reducing function first took an entire set of data for 
a given W/b and then subdivided that into sets of different 
d/a values. For each set of like d/a’s a second order least 
Mean square reduction was done using the T/a ratios as the 
input parameter. The resulting coefficients were placed in 
column form. This procedure was repeated for each different 


d/a set. The coefficients were then arranged into data sets 
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treating each location in the column as a data point. Another 
least mean square fit was done using the d/a ratios as the 
input parameter. Another set of coefficients were generated 
from this reduction. These coefficients were again placed in 
a known column order. This procedure was repeated for each 
W/b set. The coefficients were again placed in a known column 
order and the data reduced to coefficients using a least mean 
square fit with W/b as the input parameter. This procedure 
resulted in 162 coefficients. The element values could then 
be regenerated by reversing this process. The programs that 
implemented the above algorithm were tested by picking an 
arbitrary set of coefficients, generating element values then 
reducing the element values to coefficients. The coefficient 
values were recovered with differences on the order of 10°. 
This test, of course, did not give any indication of the 
performance of the algorithm with unsmoothed element values. 
The function listings for the above routines are included in 
Appendix C. 

Application of the procedure described above to the 
element values was unsatisfactory. Testing the model with 
elements generated from the curve fitting routine resulted in 
large errors. Typical error levels were found to be in the 
neighborhood of 80%. Apparently, small errors generated in 
the initial recursions were magnified by using curve fitting 
on the coefficient in the subsequent curve fits. After the 
initial attempt was found to be inadequate, several additional 
procedures were used to attempt to develop’ suitable 
coefficients. 

The element values were examined and values that were far 
from the apparent trend of the data were eliminated to create 
a data set that could be used to generate usable coefficients. 
This procedure did improve the performance but the results 


were again found to be unsatisfactory. Several possible 
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variations in the ‘correct’ trend were explored with similar 
unsatisfactory results. 

The most promising of the coefficient sets generated above 
was used to generate all of the values with the exception of 
the value of the inductor. The inductor was selected because 
it was the single most significant element to the model’s fit. 
Another program was written that found the value of the 
inductor needed, given that the other elements were all 
generated by coefficients. The elements were again reduced 
using the initial procedure described. A comparison between 
desired element values and those generated by the coefficients 
was made and showed that the elements other than the inductor 
were a perfect fit, as expected. The error in the value of 
the inductor was found to be in the 5-7% range, with some much 
larger values. This level of error in the inductor was found 
to result in 20-30% error when compared to the STRIP data. 
This procedure was applied successively to each element in the 
hope that the coefficient values would converge to an 
acceptable solution. Convergence did not occur. The output 
of the data was continually found to generate coefficients 
that produced error levels in the 20-30% range. Appendix E 
contains a portion of the MATLAB output, with typical output 
values. 

Two final attempts at function minimization were made. 
Both attempts iterated through the entire data set to allow a 
best guess at the correct coefficients. The first attempt 
used the coefficients found using the curve fitting procedure 
to generate all of the element values except the value of the 
inductor. The 27 coefficients necessary to specify the 
inductor were passed to the function minimization routine. 
This program tested every fifth data point of the entire STRIP 
data set. The error at each point was defined in the same way 
as described in the previous chapter. The total error was 
defined as the summation of the error at all the test points. 
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The resulting coefficients produced no significant change in 
the performance of the model. A modification to the program 
was made that allowed the function minimization routine to 
vary all 162 of the coefficients. This modification did not 
result in improved performance. The failure of these last two 
attempts was expected for two reasons. First, the function 
minimization routine was designed to operate with less than 
five parameters. Second, the nature of the model is such 
that when drastically incorrect element values are applied to 
the model, the value of the magnitude of, S,,, goes to either 
0 or 1 and the angle to +180°, as a result, no indication of 
increased or decreased error is given to the minimization 
routine. 

Coefficients were not found which allowed the element 
values to be represented as series of equations. Possible 
alternative approaches are discussed in the following chapter. 
The reason that least mean square curve fitting was not 
applicable is not completely clear, but it is most likely a 
combination of the following factors. 


°* The output of the STRIP program changed after the 
variables were re-specified as double precision. 
Indicating a possible problem in the program, although the 
program has been confirmed experimentally for several 
geometries. 


¢ The STRIP data may contain unlocated irregularities that 
prevent the production of smooth element values. 


°* Function minimization accentuates inconsistencies in the 
data by finding the optimum fit point by point. 


¢ Using an element generation method that results from three 
iterations of curve fitting is susceptible to significant 
error magnification during the calculation of the 
elements. 
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V. SUMMARY 


A. CONCLUSIONS 

The equivalent circuit model presented has been shown to 
accurately reproduce the response predicted by the STRIP 
program for individual geometries. Within the model, there 
are three effects that act together to replicate the behavior 
of the actual inductive strip. They are the decay of the 
field along the length of the inductive strip, the scattering 
from the edge of the strip, and the coupling of the field 
between the finline and the half-waveguide sections. The 
length of the inductive strip was accounted for by the two 
below-cutoff waveguides, which attenuated the energy flow 
along their lengths. The scattering caused by the edge of the 
inductive strip was modeled with the fixed parallel inductor 
and capacitor circuit. The coupling of the field from the 
finline to the below-cutoff waveguides was modeled with a 
frequency dependent conductance defined at the beginning of 
the strip. This frequency dependant conductance accounted for 
the transformation of the finline conductance to the 
conductance of the below-cutoff waveguides. Viewing the model 
in this way, gives valuable insight into the behavior of the 


actual inductive strip in finline. 


B. RECOMMENDATIONS 

An algorithm which accurately reproduced the element 
values was not found. To be useful in the CAD environment, a 
method of generating the element values from the geometry is 
essential. One possibility would be to smooth the STRIP data 
using some curve fitting technique prior to applying the model 
to the data. Smoothing the data first would reduce the 
tendency of the function minimization routines to produce 


element values that are difficult to describe analytically. 
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Element values that can be easily and accurately reduced to 
equations are essential because of the three levels of 
generation required. 

The model should be extended to other geometries. It has 
been tested at values of W/b as low as 0.1 and has continued 
to perform well for specific geometries. The model may be 
scalable to b/a’s of 0.4 to 0.5 as was the case for 
homogeneous’ finline, but this will require’ further 
investigation. 

The model for the loaded waveguide could be tested 
experimentally and implemented directly in CAD software. 
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APPENDIX A. DERIVATION OF LOADED HALF WAVEGUIDE MODEL 


This appendix presents a derivation of the attenuation 
through a partially loaded below cutoff waveguide and the 
computation of the equivalent height by matching the voltage- 
power impedance. The variable ‘a’ will be used here to denote 
the width of the shield. The presentation here follows 
Pozar’s development for the loaded waveguide. [Ref 7:pp. 151- 
153) 


Loaded Waveguide 





Figure 14. Loaded half- 
waveguide. 


As can be seen in Figure 14, the geometry is uniform in the y 
direction and since, n=0, the TE, modes have no y dependance 


the wave equation for H, can be written separately as 


($+ké}H,=0 (1) 





and 


(+i}H,-0 (2) 


where k, and ky are the wave numbers in the two regions (a in 
air and d in the dielectric). These two equations are related 
by the fact that both waves must propagate down the waveguide 
together. Therefore, the propagation constant, Y, must be the 
same on each side. The model assumes no loss, therefore, the 
propagation constant reduces to jf. Since, the model also 
assumes below cutoff operation B must be a positive imaginary 
number, which will result in a negative real propagation 
constant and evanescent electromagnetic waves. Therefore, B 


can be written as 
B=/e,ke-ka=/ke-k2. (3) 


Given the wave equations above and a knowledge of the boundary 
conditions, the general solution can be written as 


H, g=Acos (kyx) +Bsin (Kgx) for x<d (4) 
H,,=Ccos (k, (a-x) ) +Dsin(k, (a-x) ) ford< x<a. 


The field components of a TE field can be found from H, 





_jwop GH, _ ip Gt 
apres x,-—JP ae . (5) 


The electric field in the y direction can be found to be 


Eq 2ph (Bcos (ky x) -Asin(ky x) ) (6) 


and 
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z,.7 22 (Deos (k, (x-d) ) -Csin(k,(d-x)). (7) 


a 


Boundary conditions on the tangential electric field imply 


that B=D=0. The field components reduce to 


H,g=Acos (ky x) H,,4=Ccos (k, (a-x) ) (8) 
Ba 29RA sintkg x) Eygr TSEC sin (k, (a-2) ) (9) 


Hyg “9EA sin (kx) Hyg= ~2PE sin (ky (d-x)). (10) 


Since, the electric field tangential to the dielectric 


boundary must be continuous, the following condition applies 
E,q(d) =E,,(d). (11) 
Similarly, the tangential magnetic field must be continuous 
across the boundary 
Ha,(d) =H,,(d) . (12) 
Equations (12) and (13) can be solved and result in the 
following transcendental equation 


O=k,tan (kd) +kgtan(k,(a-d)). (13) 


Note also that 


cos (k,d) 
pA ee 4 
nA cos (Kk, (a-d)) wa 


Solving equation (3) for k, and k, then substituting into 
equation (14) produces 
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O=/ke-B? tan (dje,ke-B? ) +Ve,k¢-B? tan((a-d) /kj-B? ). (5) 


Solutions for the propagation constant can then be found 
numerically. For a homogeneous waveguide operating in the TE,, 
mode 6 can be found to be 


{eA Z) i 


Equation (15) will produce a set of f’s for different 
frequencies and geometries. By selecting a, and €,, in 
equation (18) correctly, a homogeneous waveguide can be found 
that approximates B’s frequency variation for a specific 
geometry. 


The voltage-power impedance, Z,,, is defined as 
2 
Z oe ee (17) 
Since, this is a below-cutoff waveguide, no real power will 
propagate. The imaginary voltage-power impedance can still be 


defined by using the complex power [Ref. 8]. The complex 


power can be written as 


bid a 
P= j [ Eyalxadx+ f EygHxadx |dy . (18) 
o\o d 





There is no y variation, so the integration in y can be done 


immediately 


pal feria as (19) 
0 d 
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Substituting in the known quantities 


d., - q . 
p-bf 2 Awpsin(kgx) jABsin (kgx) Ase 
0 


Kg kg 
bf -jawucos (kgd) sin (ka(x-a)) JABcos (kg) sin(k,(x-a)) 
I k,cos (ka(a-d) ) k,cos (k, (a-d) ) 
d 
(20) 
Completing the integration gives 
po A*wBpcos? (kd) tan(k,(a-d) | (a-d) A2wByp 
2k; 2kgcos (k, (a-d) ) 
a a a (21) 
_ A?oPpsin (kd) cos (kgd) , A2dw Bp 
2k 2kg 
Equation (21) can be more compactly written as 
P=bC, (22) 


where C, is the value of the terms in parenthesis in equation 
(21). Taking into account that the power is imaginary the 
average power can be written as 


=> IM(P) ; (23) 
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The voltage is given by 
b 
v= -[E, dl=-bE, (24) 
i?) 
where E, is the max field given found by finding the zero of 


the derivative of the y component of the electric field with 
respect to x 


pe 5 ee -jAwpcos (kyd) sin (k, (x-a) ) (25) 
dx\7* k,cos (k, (a-d) ) 
This results in 
O=cos(k,(x-a)) = x=(n+t) Zea ne...-1,0,1.. (26) 
2° °K 
which gives 
_ ~jAwpcos (kd) (27) 


En” kcos(k,(a-d)) ° 


The voltage-power impedance, Z,,, can then be written as 





v2 2 1 
Ve =bE, | —~ |. (28) 
oY 2P i C. 


The voltage-power impedance for the equivalent below cutoff 
waveguide can then be found as [Ref. 8]: 


2bg 4 


Zov=1207 —$——————_ . 
Aag/E-g fy (29) 
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Solving equation (29) for bs and substituting in the voltage- 
power impedance found in equation (28), defines the equivalent 
height required to produce the same voltage-power impedances 


for bot - -eguides 


_ DEny4aaf*e,q-C? (30) 


Cc. 480nf 


by 
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APPENDIX B. MATLAB FUNCTIONS FOR HALF WAVEGUIDE MODEL 


This appendix presents an explanation and listing of the MATLAB 
functions used to determine the model for the half waveguide. 
Variables defined as global_‘/name’ are global variables which must 
be specified in the MATLAB shell. The purpose of these functions 
is to produce a set of equivalent a,’s and €,,’s that approximate 
the variation of B for a specific geometry. 

The cutoff frequency was found using FIND FC.M. FIND FC.M uses 
the MATLAB function minimization routine to find the zeros of the 
transcendental equation contained in TRANS FC.M. 

The equivalent a,’s and €,4's were found using MODEL EQ.M. This 
function found the set of B’s over the frequency of interest for a 
particular geometry by finding the zero’s of the transcendental 
function found in Appendix A (equation 15). Equation (15) was 
implemented in TRANS.M. The equivalent dimensions were found by 
matching the B’s just produced with B’s calculated using test 
dimensions supplied by the function minimization routine. The 
routine selected the best choice of a,’s and €£,,'s to fit the data. 
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The following functions relate to cutoff frequency calculations: 


function fc=find_fc() 

% Finds the cutoff frequency for the loaded half waveguide. 

% The variable da refers to d/a and fc is the cutoff frequency. 
% The following variables must be declared global. global da 


tol=1.0e-5; 

for da=.01:.01:.25 
global_da=da; 
freq=13e9; % start of search 


sets the accuracy of the solution 
iterates through da 


ge de oP 


sets to global, pass to trans fc 


freq=fzero (‘trans fc’,freq,tol)% finds solution to trans fe 
fc=(fc; [ da freq/1e9]]% stores the solution 
end 


function y=trans_fc (x) 
% finds the cutoff frequency. 


beta=0; % condition for cutoff 
a=.02286/2; % wr90 in meters 
d=global_da*a; % d dimensions 
er=2.22; 

ko=2*pi*x/3e8; % wave numbers 


kd=(er*ko*2+beta%2) “*0.5; 
ka= (ko*2+beta%*2)°0.5; 

% transcendental equation -> 0 
global_yy=ka*tan (kd*d) +kd*tan (ka* (a—d) ) ; 
y=global_yy; 
end 
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The following functions relate to the computation of the equivalent 
dimensions: 


function y=model_eq(pass_ tol) 


% global variables that must be declared: global_pass tol, 
% global_aeq, global_er, global_f, global_y error global_yy 
%* g cost 


global _pass tol=pass tol; % sets tolerance throughout routine 
tol=pass_ tol; 
freq=12; % upper bound of matching 
ii=1l; x(ii,:)={.01111 1.0701];% starting point for search 
for i=.01:.01:.25 % iterates through geometries 
n=1; 
hhold=245; % beginning point for beta 
global _d=i; 
for £=6:0.5:freq % iterates through frequency 
global _f=f*le9; % pass frequency to global 
% get beta at that freq 
hhold=fzero(’trans’ ,hhold, 1le-6) ; 
g_cost (n,1)=hhold;% save the data point 
n=n+1; 
end 

% find the aeq and er that match 
x(ii,:)=fmins(’cost_eq’,x(ii,:),tol)’; 

% the beta’s held in the cost 

% matrix below. 

% save data and error 
y2(ii,:)=[i x(ii,:) global_yy global_y error ]; 
Li=iitl; 
x(ii,:)=x(ii-1,:); 


end 
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function y=trans (x) 


% This function finds the degree to which the transcendental 


% equation goes to 0 for a given beta and dimension. 


beta=x; % the value being found 

%* constants 
a=.02286/2; d=global_d*a; er=2.22; 
ko=2*pi*global_f/3e8;% wave numbers 
kd= (er*ko*%2+beta%*2)*0.5; 
ka= (ko*2+beta%*2)*0.5; 

%® transcendental equation 
global_y=ka*tan (kd*d) +kd*tan (ka* (a-d) ) ; 
y=global_y; 
end 


function y=cost_eq(x) 


global_aeq=x (1); % variables 
global _er=x (2); 


tol=global_ pass tol; % set tolerances 


freq=12; % set upper limit on matching 
hhold=245; % starting value for search 
n=1; 


for f=6:0.5: freq 
global_f=f*1e9; % passes frequency to trans eq.m 
% finds the beta that solves for 
% a particular aeg and er 
hhold=fzero(‘’trans_eq’,hhold,tol); 
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g_cost (n,2)=hhold;% saves data in matrix 


n=n+t1; 
end 
y1=0; 
for i=l:n-1 
compares the beta vectors from 


load waveguide with homogeneous 


de de de 


waveguide. error is the square 
% of the sum of the differences. 
yi=(g_cost (i,1)-g_cost (i, 2)) “2+y1l; 
end 
y_error=yl; 


y=yl; 


function y=trans_eq(x) 
% function returns the beta for a specific homogeneous geometry 
% at a specific frequency 


beta=x; % desired variable 
ko=2*pi*global_f/3e8 % calculation of wave numbers 
ka=(global_er*ko“2+beta%2) *0.5; 

yy=tan (ka* (global_aeq)); % goes to zero to satisfy B.C. 
y=yy; 

end 
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APPENDIX C. EESOF MODEL AND MATLAB MODEL 


This appendix describes the steps taken to develop the 
circuit model. The first step in the process was writing a 
MATLAB program that could replicate the performance of the CAD 
software. The CAD software available, and the target of this 
work, was the Touchstone software by EESOF. Following this 
discussion is a printout of the structure described in EESOF 
code, a typical EESOF data file and the MATLAB functions 
mentioned. 

The approach taken to develop the MATLAB code was to 
represent each of the half waveguides as a n-equivalent 
circuit, as shown in Figure 15 (the dotted lines represent 
each half waveguide element). The nature of the structure 
lends itself most easily to use conductance for calculations. 
The value of the individual elements of the t-equivalent 
circuit for a waveguide can be found as 


Y, 
= = 
Vie sinh (y,¢) Y,.=Y, tanh (y,t/2) ag 
= Ya = 
Ya" eink) Yog=Yqg tanh (ygt/2) 


where the subscript ‘d’ and ’a’ refer to the loaded and 
unloaded half-waveguides. The initial implementation of the 
above algorithm was done in a function called STRUCTURE.M. 
This function and the same model design implemented in 


Touchstone was tested over a wide range of parameter values to 


ensure that the two programs produced the same results. In 
the interest of speed STRUCTURE.M was divided into two 
functions STRUCTURE_1.M and STRUCTURE _2.M. STRUCTURE 2.M 


computed values that were dependant only on geometry and 
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n-Equivalent Circuit Representation 





Figure 15. The m-equivalent circuit for the two half 
waveguides. The dotted lines separate each of the 
waveguides. 


therefore, could be saved while the values in STRUCTURE 2.M 
depended on the element values and needed to be recalculated 
for each iteration. 

The ’structure’ functions were then called by a series of 
function minimization routines. Initially FIRST _SEARCH.M was 
used to get a rough idea of the progression of the element 
values. This program first searched over the coefficients of 
the, Y,, polynomial to get a solution reasonably close to the 
values from STRIP. Then, another function minimization was 
done using all of the parameters. The purpose of this search 


regime was two fold: a search over six parameter is very time 


41 














consuming and the results from one geometry to the next can be 
very discontinuous. By holding the value of the inductor and 
the capacitor fixed for the first search, the entire program 
ran faster and produced smoother results. After returning the 
‘best’ guess solution of the parameters to the calling 
function, the parameter values were then stored and passed to 
the next iteration as the starting point of the search. The 
function STORE.M opened a file named STORAGE.MAT and save the 
results in a column form. Figure 16 contains a flow graph 
which summarizes the above discussion. 

The performance of a function minimization is determined 
by the way the error is defined. The error, in this case, was 


FIRST SEARCH [QRRiSTRUCTURE_1 
OR 
SECOND SEARCH [SRRISTRUCTURE 2 


pom Oho k 


Oo TORE 
ELEMENTS 


Figure 16. The iteration between the functions to produce 
element values can be seen above. 











calculated, in the function CALC _ERROR.M using the square of 
the sum of the differences of the calculated values of the 
scattering parameters and the values generated by the STRIP 
program. Because both the phase and magnitude of, S,,, were 
considered equally important, the square of the difference of 
the two values was also taken and added to the error to force 
the errors to become evenly balanced. The following is an 


algebraic description error 


e,,= ( 13... -|Sy1¢l) ( ISia0}- Sarl) © (2) 


On= (LG 45-2 9430) (£9,157 29330) n (3) 


where, S,,,, is a vector that contains values of the spectral 
domain data over the frequency set and, 8,,,, is a vector that 
contains the values calculated by the model. The errors were 
combined in the following way 


Cota =Emt ept10 (@,,-p) : (4) 


After several runs of FIRST_SEARCH the values of elements 
were smooth enough to allow the calculation of equation to 
compute the element vales at specific geometries. The function 
MOD_ELEMENTS .M allowed the element values to be hand tuned to 
Manually smooth the data prior to attempting to reduce the 
data. The coefficients were calculated using FIND _COEF.M and 
FIND_WB.M. These functions were placed into another function 
called MAKE COEF.M This program iteratively called up the 
saved values of the elements and computed coefficients using 
a least mean square curve fitting routine. For a given, d/a, 
and, W/b, the second order polynomial for the, T/a, variation 
was found for each of the elements. These coefficients were 
then placed in column form. The set of data for the next, 
d/a, geometry was reduced until all the, d/a’s, for a 


specific, W/b, geometry was computed. These coefficients were 





then treated as elements and converted into polynomials that 
depended on, d/a. These coefficients are placed in column 
form and the data for each, W/b, was then reduced. These 
coefficients were again treated as elements and polynomials 
were formed using the, W/b, dependance. By reversing the 
process in the functions GET_WB.M and GET_COEF.M, the values 
of the elements were regenerated. This program was tested 
extensively by starting with an arbitrary set of coefficients 
and generating elements then re-created the coefficients and 
verifying that the original coefficients were regenerated. 

After producing a set of coefficients, the STRIP data was 
reapplied to the function SECOND SEARCH.M. SECOND _SEARCH.M 
did the same thing as FIRST _SEARCH.M but, the initial guesses 
were produced by the coefficients and the function only went 
through one minimization routine that varied all of the 
element values. Figure 17 depicts the flow through the entire 
process. 
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FIRST SEARCH [quem STRUCTURE 


ELEMENT VALUES 


COEFFICIENTS 


SECOND SEARCH 


ELEMENT VALUES 


COEFFICIENTS 





Figure 17. The entire process of developing the coefficient 
from the raw data is depicted above. 











EESOF MODEL 


! USER: Michael Linzey 


! DATE: 3 MAR 91 
! CIRCUIT: model 


! COMMENT: This is a model that represents the inductive strip 


! in finline. 


The element values must be entered for a 


! specific geometry. 


DIM 
FREQ GHZ 
RES OH 
IND NH 
CAP PF 
LNG MIL 
TIME PS 
COND /OH 
ANG DEG 


VAR 

! GEOMETRY SPECIFIC 
a=900 
b=400 
Ww=200 
d=0 
T=10 

! ELEMENTS 
1_d=12.6320 
c_d=.0064 
turn0=1.6020 
turn1=3.1142 
turn2=-2 .5861 
turn3=0.1429 


EQN 
PI=3.14159 
A_a=a/2 


} MODEL 
balance=1 
ba=b/a 
wb=W/b 
dal=d*2/a 
da=d/a 


! EQUIVELENT VALUES FOR THE LOAD HALF WAVEGUIDE 
er dl=.999851+dal* .0149217-dal1**2* .3391197 
er _d=er dl+dal**3*10.61754-dal**4*3.664746+da1l**5*.135657 
A_dl=1+dal* .001314-dal**2* .0302559+dal**3* .342420 
A_d=(A_dlil-dal**4*2.161754+dal**5*1.7460658) *a/2 
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c=2.998*10**8*1000/ .0254 
mu=4*pi*10** (—7) /(1000/.0254) 
£=FREQ*10**9 
k0=2*pi*f/c 

! beta presumed imaginary 
beta=sqrt (abs (er_d*k0**2-(pi/A_d) **2) ) 
kd=sqrt (2.22*k0**2+beta**2) 
ka=sqrt (k0**2+beta**2) 


CS=cos (d*kd) 

SN=sin (d*kd) 

zov_n=2*b*ka*kd**3* (2*pi*f) *mu*CS**2 

zovl=ka**3* (SN*CS-d*kd) *cos (A _ d*ka-d*ka) **2 
zZOov2=kd**3*CS**2* sin (A _ d*ka-d*ka) *cos (A_ d*ka-d*ka) 
ZOV _ d= (beta* (zovl+zov2-ka*kd**3* (A _ d-d) *CS**2) ) 
zov=zov_n/zov_d 


B_d=-zov*sqrt (abs (4*A_d**2*er d*£**2-cx*2) ) / (480*f*pi) 
! CALCULATES THE FINLINE EQUIVALENT DIMENSIONS 
c1=-4.9723*ba**2+4.7413*ba-0.7651 
Aeg= (2- (1- (bat .45) * (l-wb) **2) **0.5+c1* (l-wb) **26) *a 


c9=-20.16*da**2+6.42*dat+.6494 
erl=c9* (1-wb** (l-exp (-10*da) ) ) 
er=erl+wb+2 .604*da+ (1—da) **6* (1-wb) 


c2=(-115.79*da**2+27.87*da-—. 4933) *ba+87.52*da**2—-22 .49*da-.1932 
c3=0.29+0.0773*exp (1-40*da) 
c4=(20.1154*da**2-3.729*da-0 .0611) *bat (~26.1788*da**2+5 .537*datl1.0376) 
c5=-13.5217*da**2+2.4017*da+0.0411 

Beq avel=c2* (1l-wb** (2*ba*c3) ) +c4+c5* (l-abs (ba-wb) **2) **4 
Beq_ave=(Beq avel-0.025* (1-abs (.925-—wb) **2) **16) *b 

c6= (-76.251* (da) **24+17.23*da-.1578) *bat111.2*da**2-20 .84*da-—.2936 
c7=(64.82*da**2-14.77*da—.3029) *ba-107 .1*da**2+22 .85*da-.2936 
c8=(9.696*da**2-1 .449%da—-.1431) *ba-12.13*da**24+1.39*da+.1195 

m=c 6 *wbh* *2+c7 *wb+c8 

fco=(29980/2.54) / (2%a) 

Beq=m* (FREQ/fc~1.56) *b+Beq_ave 


! IMPEDANCE TRANSFORMER 
Lambda= (29980/2.54) /FREQ 
LpovL1=1/ (1- (Lambda/ (2*Aeq*er**0.5) ) **2) **0.5 
LpovJ.2=1/ (1- (Lambda/ (2*A) ) **2) **0.5 
Z1=i120*PI* (2*Beq/Aeq) *LpovL1l/er**0.5 
f n=(FREQ-fc) /fc 
Zz _ d=1000/abs (turn0+turnl*f | ntturn2*f n**2+turn3*f_n**3) 
Zs  a=Z __d/balance 


X_a=(Z1/2_a)**0.5 
X_d=(Z1/Z_d)**0.5 
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CKT 


IND 10 L*Ld 

CAP 10 c*’CcC d 

RWG 12 A*A_a B“*B 1“°T ER=1 RHO=1 
IND 20 L*L d 

CAP 20c*%c d 

DEF2P 1 2 Air 

IND 10 L*Ld 

CAP 10c*cd 

RWG 1 2 A*A_d B*B d L*T ER*er_d RHO=1 
IND 20 L*‘Ld 

CAP 20 c*c_d 

DEF 2P 1 2 Dielec 

XFER 12 00 N*X_a 

XFER 13 00N*%X a 

Air 24 

Dielec 3 5 

XFER 6500N*X d 

XFER 6 40 0 N*X_a 

DEF2P 1 6 STRIP 

RWGT 1 A*Aeq B“Beq ER“er RHO=1 
DEF1P 1 WEDGE 


!SPECTRAL DOMAIN DATA 
S2PA 120 t50d30w5.S2P 
DEF2P 1 2 SPEC_DOM 


SPEC_DOM 1 2 
DEF2P 12 SPEC_M 


TERM 
STRIP WEDGE WEDGE 
PROC 
OUT 
STRIP MAG[{S11] GR1 
STRIP ANG[S11]) GRIA 
SPEC _DOM MAG[S11) GR1 
SPEC DOM ANG[S11] GRI1A 
FREQ 
SWEEP 8 12 1 
GRID 
RANGE 8 12 1 
GR1 0 1 -1 
GRIA 80 180 


48 








The following listing illustrates the form of the spectral domain input 
data file for the iductive strip model. This is the listing called by 
the model on the previous page. 


!Spectral Domain data file for the following geometry 
{ T=50 d=30 wW=200 


'Freq IS11| <S11 }$12| <S12 }S21] <S21 1S22|] <S22 
8.0 0.8992 146.6191 0.1914 56.6191 0.1914 56.6191 0.8992 146.6191 
9.0 0.8528 139.1836 0.2727 49.1836 0.2727 49.1836 0.8528 139.1836 
10.0 0.8056 132.0117 0.3510 42.0117 0.3510 42.0117 0.8056 132.0117 
11.0 0.7568 126.7910 0.4272 36.7910 0.4272 36.7910 0.7568 126.7910 
12.0 0.7093 121.3066 0.4968 31.3066 0.4968 31.3066 0.7093 121.3066 
13.0 0.6610 115.7168 0.5630 25.7168 0.5630 25.7168 0.6610 115.7168 
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The following is a listing of MATLAB functions used in element calculations: 


function [(f,mag,ang]=structure(fl,a,b,d,W,t,1 d,c_d,t0,t1,t2,t3) 
% This function replicates the performance of TSTONE software for the 
% configuration under consideration. 
% calculates constants 
a_a=a/2; er_a=1; wb=W/b; 
fc=2.998e8*1le-9*1000/ .0254/ (2*a) 
% iterates through frequencies 
for n=1: length (f1) 
freq=f1(n) ; % get first freq 
f=freq*1e9; 
% get equivalent loaded waveguide 
{a_d,b d,er_dJ=half_guide(a,b,d,f); 


y0O_a=yov(a_a,b,er a,f); % calculate the susceptance of 

yo _d=yov (a: a,b _ d,er d,f); % the two guides 

gamma | a=real (j*betatf, a_a,er_a)); % calculate the prob. const. 
gamma_d=real (j*beta(f, a_d, er a)); 


fn= (freq-fc) /fc; % calculate the susptance used in 
yO_f=0.001*abs (t0+t1l*fn+t2*fin*2+t3*fn%*3) ;% turns ratio 


yO_a=y0_ a/y0 f; % normalizes the impedance 
yO d=y0 d/y0 f; 

% pi equivilant for waveguide 
ya=y0_a/sinh (gamma_a*t) +y0_d/sinh (gamma_d*t) ; 
yb=y0_a*tanh (gamma_a*t/2) +y0 _d*tanh (gamma_d*t/2) ; 

% tank cuircuits 
tank_a=(j*2*pi*f*c_d+1/(j*2*pi*f£*l_d))/y0 £; 
tank_d=(j*2*pi*f*c_d+1/(j*2*pi*f*l_d))/y0_ £; 


yb=yb+tank_| at+ttank _d; % form the complete pi =quivilant 
y_eq=1/(1/(1+yb)+1/ya)+yb; % find the equivilent suseptance 
sll= (1l-y_eq) / (1l+y_eq) ; % calculate sli 


g(n,:)=[f*le-9 abs(sll) angle(sil) *180/pi); 


end 
f=g(:,1); 
mag=g (:,2); 


ang=g(:,3); 
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function structure _1(f1,a,b,d,W,t) 

This function coupled with STRUCTURE 2 produce the same results as 
STRUCTURE. STRUCTURE_1 and STRUCTURE 2 together repoduce the 
performance of TSTONE for the device under concideration. STRUCTURE_1 
computes the geometry specific parts of STRUCTURE to prevent 

for each frequency iteration. 


de de dP dP de 


% constants 
a_a=a/2; n=0; er_a=1; wb=W/b; 
fc=1le—-9* (2.998e8*1000/.0254) / (2%a); 

% iterate through the frequencies 
for ifreq=1:length (f1) 


freq=f1 (ifreq) ; % select first freq (in GHz) 
f=freq*1le9; 
{a_d,b_d,er_dj=half guide(a,b,d,f); % get the equivilent loaded guide 


yO a (ifreg) =yov(a_ a, ib; er a,f); 

yO dad (ifreq) =yov (a: d,b a, ex! d,£); 
gamma_ a=real (j*beta(f,a_ a,er a)); 
gamma_d= real (j*beta(f,a | d,er a) 3 
gamma_check (ifreq)=gamma_d; 


calculate the impedances for 
the two guides 
calculate the probigation 
constants 
used to check if guide starts 
% to probigate 
sinh __a(ifreq) =sinh (gamma _| a*t); % calculates need trig fucntions 
sinh | _d(ifreq) =sinh (gamma _ a*t); 
tanh_a(ifreq) =tanh (gamma_a*t/2) ; 
tanh | d(ifreq) =tanh (gamma_ ~da*t/2); 
end 
f=f1; 


dP dP dP dP dP 


% saves calculations to CONSTANTS 
save constants f yO a yO d gamma_check fc sinh_a sinh_d tanh_a tanh _d 
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function [f, mag, ang])=structure_2 (1 _d,c_d,t0,t1,t2, t3) 
% This function completes the calcuations started in STRUCTURE_1; includes 
% the effects of the elements. 
load constants % gets constants from STRUCTURE 1 
1_d=abs(1_d); c_d=abs(c_d); n=1; fl=f; % takes the abs of L and C to 
% prevent the use of negative #’s 
for freqi=1: length (f) % iterates through frequency 
if gamma_check (n) ==0, % checks to see if the guide is 
% probigating, skips if yes 
else 
freq=fl(freqi); f=freq*le9; w=2*pi*f; % variables 
% calculates transformed impedance 
f n=(freq-fc)/fc; 
yO_f=0.001*abs (tO+tl*f n+t2*f n*2+t3*f n%3) ; 
% Normalizes the suseptance of 
% guides 
yO_a(n)=y0_a(n)/y0 f; 
yo. ~d(n)=y0 ~d(n)/y0_ f; 
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% pi-equivelant of guides 
ya=y0_ a(n) /sinh_a(n)+y0_d(n)/sinh_d(n); 
yb=y0_ a(n) *tanh a(n) +y0_ a(n) *tanh © . a(n); 

% suseptance of the tank circuits 
tank a=(j*w*c_ d+l1/(j*w*l_d))/y0_ f; 
tank_d=(j*w*c  d+1/ (j*w*l_d)) /y0_£; 

% total circuit 
yb=yb+tank_attank d; 

% egqivilent suseptance 
y_eq=1/(1/ (1+yb) +1/ya) +yb; 

% calculate s1l 
sll=(1-y_ eq) / (ity _eq); 

g(n,:)=[£*le-9 abs(sll) angle(si1)*180/pi]; 
Ean 
end 
end 
=9(:,1); 
mag=g(:,2); 
ang=g(:,3); 


The following functions support the ’structure’ functions: 


function b=beta(f,a,er) 
tcalcualtes beta for a given frequency for the TE10 mode. Assumes that 
%a is measured mils. 


c=2.998e8*1000/ .0254; 
k=2*pi*f£* (er) *.5/c; 
kc=pi/a; 

b= (k*2-ke%2)%*.5; 

end 


ne re ree nc ee ee re ee ee ee ee ee ee ee ee 0 ee re ee rr ee ee re ee ne ee ee ee ee ee ee ee ee ee oe eee ee 


function {au d,b_d,er_d]=half guide (a,b,d, f) 
%This function computes the effective dimentions of a the below cutoff 
twaveguide section that contains dielectric. The model assumes that 
%the dielectric is er=2.22. 
% constants 
c=2.998e8*1000/.0254; mu=pi*4e-7/ (1000/.0254); d=d*2/a; 
% calc er 
er d=.999851+d*.0149217-d*2* . 3391197 +d*3*10. 617549-d*4*3.664746+da*5*.135657 


a_d=a/2* (1+d* .001314-d*2* .0302559+d*3* .342420-—d*4*2 .161754+d°5*1.7460658) ; 
% calculate Zov 

kO=2*pi*f/c; 

beta=sqrt (er d*k0*2-(pi/a_d)%2); 

kd=sqrt (2. 22*k0*2-beta*2) ; 

ka=sqrt (k0*%2-beta%2) ; 
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CS=cos (d*kd) ; 

SN=sin (a*kd) ; 

zov_n=2*b*ka*kd*3* (2*pi*f) *mu*CS*2; 

zovi=ka*3* (SN*CS—d*kd) *cos (a_d*ka-d*ka) *2; 
zov2=kd*3*CS*2*sin (a_d*ka-d*ka) *cos (a_d*ka-d*ka) ; 
zov_d=(beta* (zovl+zov2—-ka*kd%3* (a_d-d) *CS*2) ); 
zov=zov_n/zov_d; 

b_d=real (-zov*sqrt (4*a_d*2*er_ d*£*%2-c%2) / (480*£*pi) ) ; 


a a es ee ce ae re es ee ee ce cee ce ee ee ee en ee ee ee ee ee ee ee ee ee ee ee ee er ee ee ee ee 


function y=yov(a,b,er, f) 

% calculates voltage power suseptance 
fc=(1000/ .0254) *2.998e8/ (2*a*er*0.5) ; 
lambda= (2 .998e8*1000/ .0254) /f; 
z=(376.7/er*0.5) *(2*b/a) / (1— (f£c/£) *2) *0.5; 
y=(1/z); 


53 











The following functions are connected with the search routines: 


function error=calc_error(mag,s_mag,phase,s phase, f, x) 
% This function caclulates the error between the two input magnitudes 
% and phases. 
s_Mag=s_mag (1: length (mag) ) ; % sets the length of the vectors to be 
s phase=s_ phase (1:length (mag) ) ; % the same, needed in the case of 
% propagation less than 13GHz 
m= (mag-s_mag) ; 
p=(phase-s_ phase) ; 
for i=1: length (m) % finds the relative error at each point 
m(i)=m(i)/s_mag(i); 
p(i) =p (i) /s_phase (i) ; 
end 
global_m=max (abs (m) ) *100; 
global _ p=max (abs (p) ) *100; 
error mag=m’ *m,; 
error ph=p’ *p; 





% calculates the error returned 
error=error mag+terror ph+10* (error _mag-error_ ph) “*2; 
% viewing options when running in 
% foreground 
%[global_m global p x’] 
splot (f,s_mag,’r’',f,s phase/180,’g’,f,mag,’*r’,f,phase/180,'*g’) 


SS A YO A TS SS CY SS SS AD SS SS AS SY RY OS ce AD SS SS SS SS SS SD SS SS A SS OS 


function error=model (x) 

This function takes input from the matlab function minimization 
routine passes values to STRUCTURE 2. STRUCTURE_2 calculates the 
response and return the values of the S11 scattering parameters. 
CALC_ERROR then compares the two responses and calculates an error 
which is returned to the function minimization routine, for use in the 
production of the next guess. 


de dP dP OP HP CP 


load sil _data_m 

1_d=x (1) *le-9; % element values 

c_d=x (2) *le-12; 

turn0=x (3); 

turnl=x (4); 

turn2=x (5) ; 

turn3=x (6); 
% finds response for given element values 

(f mag phase]=structure_2(1_d,c_d,turn0,turnl,turn2,turn3); 
% calculates error 

error=calc error (mag,s mag,phase,s phase, f,x) ; 


te cae ee cm es er eat ee seme Se cme ee ee cee ee ce ey ce ee ee ee ee ee ee ee ws ee ee ee ee ee ee ee ee ee oe 


function error=model 2 (x) 
% This function takes input from the matlab function minimization 
% routine passes values to STRUCTURE_2. STRUCTURE 2 calculates the 
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% response and return the values of the S11 scattering parameters. 
% CALC ERROR then compares the two responses and calculates an error 
% which is returned to the function minimization routine, for use in the 
% production of the next guess. 
load sll _data_m 
load g_ hold 
1_d=g_ hold(1) *1le-9; % element values 
c d=g | y hold (2) *le-12; 
turn0=x (1) ; 
turnl=x (2) ; 
turn2=x (3); 
turn3=x (4); 
% finds response for given element values 
[f mag phase]=structure 2(1 d,c_d,turn0,turnl,turn2,turn3) ; 
% calculates error 
error=calc_error(mag,sS_mag,phase,s phase, f,x); 


function g=first_search(g init,tol,begin) 
This function finds and stores element values that match the scattering 
parameters found in the file S11 DATA. The data in S11_DATA is in 
column form as follows: 

{ freq, sll_mag, Sll_phase, t, d, W] 
The data is first broken in to blocks that represent each geometry. The 
data is then saved in S1l1_data_m for recall by the function minimization 
routines. STRUCTURE 1 then calculates the geometry dependant parameters. 
Initail starting points for element values are eather given by the 
argument or from the previous iteration. Function minimiztion routines 
are then called to find values of the elements that match the S1l_DATA. 


oe dP oP dP dP dP CP OP FP oP 


The variables global_m and global_p must be declared to be global. 


%* initial values 

g=g_init; g_select=1; d_hold=0; a=900; b=400; n=0; 

load sll _data 
% block data by geometry, using the 
% starting frequency 

for i=l:length(sll data) 


if 7 == $11 _data(i,1), 
n=n+1; 
block (n) =i; % record begining of each block 
end 
end 
block (n+1)=length(sll1_ data) +1; % mark the end of the last block 
for i=begin: length (block) ~1 % iterate through all the data 
i_end=block(it+1)-1; % calculate the end of the block 
f=sll _data (block (i) +1: i end,1); % block out the data 


s _mag=s11 _data (block (i) +1: i_end, 2); 
s_phase= =s1l _data (block (i) +1: i_end, 3); 
t=sl1l _data (block (i) , 4); 

d=s1l1_ data (block (i),5); 
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W=sll_ data (block (i),6); 
i_end=length (s_mag) ; 


for k=length(s_mag)-1:-1:2 % test to see if magnitude of 
if s_mag(k)<s_mag(k+1) % sll is increasing 
i_end=k; 
end 
end 
f=f (1:i_end) ; % eliminate the bad points noted 
s_mag=s mag (1: i _end) ; % above 
s_phase=s_phase (1: i _end) ; 
{d t W] % display block being worked on 


(£ s_mag s_phase] 
save sll data m 


structure 1(f,a,b,d,W,t) ; % calculate geometry dependant 
% values 
g_hold=g % save guess of element values 


save g hold g hold 
~ ~ % pass turns ratio into function 
% minimization routine 
g_hold=fmins (’model_ 2’, [{g(3) 9(4) g(5) 9(6)}],0.10); 
~ % re-form element vector 
g=[g(1) g(2) g_hold(1) g_ hold(2) g_hold(3) g_hold(4)] 
g_hold=g = 
[global | m global p] 
save g . hold g_hold 


[g, count ]=fmins (’model’ ,g, tol) ; % pass all elements into function 
% miniminzation routine 
store (a,b,W,d,t,tol,g’ ,count/100) % save the element values 
if t<15, 
g_save=g % save the value of the begining 
end % of one d set to start the next 
if t>950 
g=g_save 
end 
end 


function g=second search (tol, begin) 

% Operation of this function is the same as the FIRST SEARCH except 
% in the way the initial guess of the elements is found. As note 

% below. 


n=0; a=900; b=400; 
load sll data 
for i=l: length(sll_data) 
if 7 == sll_data(i,1), 
n=n+1; 
block (n) =i; 
end 
end 
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block (n+1)=length(s11_data) +1; 
for i=begin: length (block) -1 
i_end=block (i+1)-1; 
=311 _data (block (i) +1: i _end,1); 
3 _mag=sil | data (block (i) #1:i _end, 2); 
s phase=sil | data (block (i) +1:i _end, 3); 
t=sll_ data (block(i),4); 
d=s11_ data (block (i), 5); 
W=311 data (block (i) ,6); 
i_end=length(s_mag) ; 
for k=length(s_mag)-1:-1:2 
if s mag(k)<s mag(k+1) 
i end=k; ae 
end 
end 
f=f (1:i_end); 
s_mag=s "mag (1: i end) ; 
s_phase=s_ phase (1: i_end); 
{d t W] 
{f s_mag s_phase] 
save sll data _m 
structure _1(f,a,b,da,W,t) ; 
% gets the elements as a function 
% of parameters 
[L,C,t0,t1,t2,t3]=get_elements(a,b,d,W,t) ; 
g=(L C tO t1 t2 t3) 
[g, count ]=fmins (’model’,g,tol) ; 
store2 (a,b,W,d,t,tol,g’,count/100) 
if t<15, 
g_save=g 
end 
if t>950 
g=g_save 
end 
end 
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The following functions relate to coefficient reduction and element 
generation. 


function coef d=find_coef (s) 
% This function find the coefficients of the polynomials using MATLAB’s 
% mean square function. 
Nt=2; Nd=2; Nb=2; 
$(10:11,:)=abs(s(10:11,:)); 
s=[s zeros (16,1)]; 
for n=1:length(s(5,:))-1 
if s(5,n)>s(5,n+1) 
Ne=n; 
for i=10:15 % finds the coef of the elements 
% with respect to T/a 
coef t=[coef_t; polyfit(s(5,Nb:Ne)/s(1,1),s(i,Nb:Ne),Nt)’ ] ; 


end 
coef_l=[(coef_1 coef_t] % puts the results in a column 
d={d | 3 (4, Nb) /7s(1,1)] % makes d/a list to reduce the 
% above list 
coef t=[]; 
Nb=Ne+1; 
end 

end 

[row, col]=size(coef_1); 

for n=l:row % finds the coef of the coef 


% generated above wrt d/a 
coef_d=(coef_d; polyfit(d,coef_1(n,:},Nd)’ ] 
end 
function coef wbh=find_whb (wb, coef_in) 
% Finds the coefficients with respect to W/b 
(xr, cJ=size(coef_in); 
Nw=2; 
for i=l:r 
coef wh=[coef_wb; polyfit (wb,coef_in(i,:),Nw)]; 
end 
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function [L,C,t0,t1,t2,t3]=get_elements (a,b,d,W,t) 
% This function used the values of the coefficients in the file COEF 
%* to calculate the element values. 


load coef 


Nd=2; % the order of the d/a polynomial 
% needed for initial decomposition 
coef _d=get_whb (Nt, coef, W/b) ; % accounts for the W/b variation 


{x,c]=size (coef d); 

N=r/6; da=d/a; 

for i=1:N % iterates throught the coef; accounts for 

% the d/a variation 

coef_L=[ coef_L polyval(coef_d(i,:),da) ]; 
coef | _C=[ coef C polyval (coef _. a(N+i,:),da) j; 
coef _t0=[ coef t0 polyval (coef | a(2*N+i, : ),da) 
coef tl=[ coef tl polyval (coef _ a(3*N+i,:),da) 
coef t2=[ coef t2 polyval (coef _ da (4*N+i,:),da) 
coef | ~t3=[ coef | ~t3 polyval (coef _ a(5*N+i,:),da) 

end 


Sa 


me Me Me Ne 


% computes the value of the elements 

=polyval (coef _L,t/a); 

C=polyval (coef | C,t/a); 

t0=polyval (coef __ to, t/a); 

tl=polyval (coef t1,t/a); 

t2=polyval (coef _- ~t2,t/a); 

t3=polyval (coef _t3,t/a) ; 

function coef_d=get_wb(Nb, coef wb, wb) 

% Account for the W/b variation 

for i=l: length (coef_wb) 
c0(i)=polyval (coef wb(i,:),wb); % calculates the W/b variation 

end 


for i=1:Nd+1: length (c0) % resizes the matrix for next 
cl=[ cl; cO(i:it+Nd) J]; % operation 
end 


coef _d=cl; 
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APPENDIX D. COEFFICIENT REDUCTION OUTPUT 


Following this page are three sample outputs’ that 
illustrate the difficulties encountered during the development 
of the coefficients which were to represent the element 
values. The first listing shows the results of a function 
minimization where the inductance was allowed to vary and the 
capacitance and the four terms of the turns ratio were 
calculated by the coefficients. The second listing shows the 
percent error in the desired value of the inductance and the 
value produced from the least mean square generating programs. 
The third listing shows the percentage error in phase and 
magnitude with the STRIP data when the coefficients calculated 


above were used to generate element values in the model. 


60 





The following listing resulted from generating element values for 


all of the elements except the inductor using coefficients and then 
using a function minimization routine to generate the inductor 


value. The columns are arranged as follows: 


HAaston 


tolerance 
(search count) *0.01 
% phase error 
% mag error 
L inductor in nH 

C capacitor in pF 
tO O order coefficient for the turns ratio 
tl 1st order coefficient for the turns ratio 
t2 2nd order coefficient for the turns ratio 
t3 3rd order coefficient for the turns ratio 


Columns 1 through 7 


_ 
WRNHRFOTIONOO 


.0000 900.0000 900.0000 900.0000 900.0000 900. 


b> bb 
a oo 
oO oo 


.0000 400.0000 400.0000 400.0000 400.0000 
.0000 400.0000 400.0000 400.0000 400.0000 

0 0 0 0 0 
0000 20.0000 30.0000 40.0000 50.0000 
0100 0.0100 0.0100 0.0100 0.0100 
0900 0.0900 0.0900 0.0900 0.0900 
6978 1.8785 0.9913 0.9183 0.8394 
4667 0.5811 0.8052 1.0324 1.2725 
6283 17.0645 16.7005 16.4661 16.3048 
0047 0.0049 0.0052 0.0054 0.0057 
5539 1.4595 1.3682 1.2798 1.1946 
5056 2.9212 3.3246 3.7159 4.0949 
5626 0.5833 -0.3683 -1.2922 —-2.1885 
2824 -2.6831 -2.1004 -1.5344 -0.9851 

0 0 0 0 0 
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ay 


OWPRFODFF OO 


Columns 8 through 14 


900.0000 900.0000 900.0000 900.0000 900.0000 900.0000 900.0000 
400.0000 490.0000 400.0000 400.0000 400.0000 400.0000 : 
400.0000 400.0000 400.0000 400.0000 400.0000 400.0000 400.0000 


P< 
oO 
i=) 
So 
oO 
o 
[~) 


0 0 0 0 0 0 0 
80.0000 90.0000 100.0000 200.0000 300.0000 400.0000 500.0000 
0.0100 0.0100 0.0100 0.0100 0.0100 0.0100 0.0100 
0.0900 0.0900 0.0900 0.1000 0.1000 0.0800 0.0800 
2.5424 1.1586 1.1907 0.6524 0.9807 3.7443 2.2694 
2.3146 0.3406 0.6124 0.9522 0.8626 0.0775 0.2237 
15.8997 16.1901 16.1876 16.5842 15.8695 13.2978 9.6489 
0.0065 0.0068 0.0071 0.0101 0.0136 0.0174 0.0216 
0.9570 0.8838 0.8137 0.2799 0.0499 0.1236 0.5011 
5.1589 5.4892 5.8073 8.3180 9.6100 9.6833 8.5378 
-4.7114 -5.4970 -6.2549 -12.3128 -15.6042 -16.1289 -13.8871 
0.5630 1.0457 1.5117 5.2560 7.3344 7.7469 6.4936 
0 0 0 0 0 0 0 


Columns 15 through 21 


900.0000 900.0000 900.0000 900.0000 900.0000 900.0000 900.0000 
400.0000 400.0000 400.0000 .0000 400.0000 400.0000 400.0000 
400.0000 400.0000 400.0000 400.0000 400.0000 400.0000 


» bd 
oo 
oo 
oO 
Oo 
Oo 
foo] 


10.0000 10.0000 10.0000 10.0000 10.0000 10.0000 10.0000 
10.0000 20.0000 30.0000 40.0000 50.0000 60.0000 70.0000 
0.0100 0.0100 0.0100 0.0100 0.0100 0.0100 0.0100 
0.0900 0.0900 0.0900 0.0900 0.0900 0.0900 0.0900 
2.7920 19233 1.2957 0.5927 0.5989 0.5894 0.7307 
1.0093 1.1187 1.2459 BR 7 1.4925 1.6206 1.7332 
15.4365 14.7504 14.2990 13.9707 13.7161 13.5044 13.3290 
0.0047 0.0048 0.0049 0.0050 0.0051 0.0052 0.0054 
1.4917 1.4271 1.3647 1.3045 1.2465 1.1907 1.1370 
3.2785 3.4819 3.6807 3.8750 4.0648 4.2501 4.4309 
=0.<6303 -1.1668 -1.6912 —-2.2036 -2.7039 3.1922 -3.6683 
~1.8272 -1.4472 -1.0769 -0.7162 =-0.. 3651 -0.0237 0.3081 


0 0 0 0 0 0 0 





Columns 22 through 


900.0000 
400.0000 
400.0000 
10.0000 

80.0000 

.0100 
.0900 
-7781 
.7554 
- 1603 
.0055 
.0855 
.6071 
.1324 
. 6303 
0 


~ 
OhPhRFPOWF OC O 


el 


Oh RHF OWH OC O 


Columns 29 through 


900.0000 
400.0000 
400.0000 
30.0000 
10.0000 
.0100 
.0900 
.0938 
-4272 
-4914 
.0049 
-4964 
. 3038 
.7174 
1571 

0 


es 
OWfFOUOrFOO 


900. 
400. 
400. 

30. 

20. 
.0100 
.0900 
.0476 
.6167 
.9184 
.0048 
.4717 
-2170 
.6376 
.2234 


wp 


OWhrRrRPOBDOrF OO 





0000 
0000 
0000 
0000 
0000 


0 


= 


FUP OONF OOO 


WwW 
oO 


900. 


& 
oO 
oO 


400. 


GW) W 
oo 


e 
OWhrRPOPROrF OO 


a 


WoOHnHODOrFWOOSO 
feo] 
jo) 
~ 
Ne) 
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f=] 


| 
ae 
OQrATIOOOrFOOO 


.0000 
.0000 
.0000 
.0000 
.0000 


i 
od 
ONAITIOAOWOFOOO 


900. 


r= 
o 
oOo 


400. 


ao WwW 
oo 


ray 


OWWrROROOOSO 


1 
- 
UNATIOONOOCOO 





Columns 36 through 42 


.0000 
.0000 
.0000 
. 0000 
.0000 
-0100 
.0900 


.0000 
-0000 
- 0000 
. 0000 
-0000 
.0100 
.0900 
-9504 
.1042 
-8818 
.0045 
. 3245 
-7701 
. 3334 
.6914 


0 


Columns 43 through 49 


.0000 
-9000 
.0000 
. 0000 
.0000 
.0100 
.0900 


0.5845 


- 9808 
-9465 
-0049 
- 6621 
4373 
~7357 
- 8258 


0 


.0000 


roy 
OP BF ORF OOO 
WwW 
fon) 
ra 
fo) 
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NAF OWODOCOOC SO 


900. 


OWWFrFOL- 





1 
rary 
WORDrFONWOOCO 


900. 


> b& 
Noo 
ooo°o 


a 


OWWrFOoORO- OO 





.0000 


0000 
0000 
0000 


.0800 
. 1880 








Columns 50 through 56 


.0000 
- 0000 
-0000 
.0000 
.0000 
.0100 
.0900 


1.1351 


rR 
OPeNOUOOOCOSO 


-7281 
-9922 
-0047 
-6004 
-4442 
.0665 
- 4303 


0 


r= 


OphroOROOCO°O 


.0000 
.9000 
.0000 
.0000 
.0000 
-0100 
.0900 


1.1735 


- 4689 
“9515 
.0047 
.5910 
. 3108 
.8475 
.3793 


0 


.0000 
.0000 
.0000 
.0000 
.0000 
.0100 
0900 
- 4383 
. 7000 
.6776 
-0053 
- 9983 
-1426 
-4611 
. 3600 


0 


= 


ObPHhRrF OO. 0000 





900.0000 
400.0000 
400.0000 
50.0000 
200.0000 
0.0100 
0.0800 
1.2297 
0.5446 
3.4752 
0.0060 
1.5040 
2.4512 
1 
0 


65 


ray 


OPWROOOCO0O° 


.0000 


0.0800 


FoMmro~TIN Oo 


900. 


» 
oO 
oOo 


400. 


~~! 
ooo°o 


ay 


COWWFROWOrF 


~9515 


Columns 64 through 70 


900.0000 900.0000 900.0000 900.0000 900.0000 900. 
400.0000 400.0000 400.0000 400.0000 400.0000 400. 
400.0000 400.0000 400.0000 400.0000 400.0000 400. 
70.0000 70.0000 70.0000 70.0000 70.0000 70. 
80.0000 90.0000 100.0000 200.0000 300.0000 400. 


0.0106 0.0100 0.0100 0.0100 0.0100 0 
0.0900 0.0900 0.0900 0.0800 0.0800 0 
1.1151 1.2713 1.2846 1.2579 0.7262 0 
0.6889 0.6666 0.6126 0.6307 1.7171 2 
13.3299 13.2318 13.1585 12.5084 11.8911 10 
0.0063 0.0065 0.0067 0.0094 0.0134 0 
1.8584 1.8357 1.8133 1.5995 1.4051 1 
3.4892 3.4009 3.3185 2.8174 2.9043 3 
-3.2833 -3.1269 -2.9819 -2.1580 —-2.4734 -3 
0.0488 0.0064 -0.0332 -0.2808 ~0.2573 0 
0 0 0 0 0 


~] 
~J 


Columns 71 through 


900.0000 900.0000 900.0000 900.0000 900.0000 900. 
400.0000 400.0000 400.0000 
400.0000 400.0000 400.0000 400.0000 

90.0000 90.0000 -0000 90.0000 


r=) 
oO 
co) 
oOo 
i=) 
Oo 
oO 
oO 
oO 
fo) 
oO 
& 
o 
Oo 


400. 


b> 
nooo 
ooo°o 
oo 
oo 
oo 
oo 
© 
oO 


Ww 
oo 
fo) 
f=) 
oO 
co) 


10.0000 20.0000 40.0000 0000 60 
0.0100 0.0100 0.0100 0.0100 0.0100 0 
0.0900 0.0900 0.0900 0.0900 0.0900 0 
2.2465 0.8248 0.6130 1.0333 1.4566 i 
1.1545 1.2252 1.3701 1.4892 1.6555 1. 

15.1502 14.2180 13.5937 13.1212 12.7360 12. 
0.0054 0.0059 0.0064 0.0070 0.0075 0. 
2.5436 2.4804 2.4179 2.3561 2.2952 24 
3.2151 3.3331 3.4468 3.5563 3.6615 co 

-2.5900 ~2.8140 -3.0304 —-3.2393 —-3.4408 =3.. 

~1.2988 ~1.1741 -1.0538 -0.9381 -0.8268 =0. 

0 0 0 0 0 
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ONOfrhRrFOMNDOSO 


ay 


OWWNONNFHOO 











Columns 78 through 84 


900.0000 900.0000 900.0000 900.0000 900.0000 900. 
400.0000 400.0000 400.0000 400.0000 400.0000 400. 
400.0000 400.0000 400.0000 400.0000 400.0000 400. 
90.0000 90.0000 90.0000 90.0000 90.0000 ; 
80.0000 90.0000 100.0000 200.0000 300.0000 400. 


Oo 
oO 


0.0100 0.0100 0.0100 0.0100 0.0100 0 
0.0800 0.0800 0.0800 0.0800 0.0800 Q. 
3.5099 3.7978 3.9986 3.0324 2.1798 0. 
1.4109 1.5650 1.7247 -0.6245 0.1249 2. 
12.5358 12.2834 12.0626 10.5341 9.7944 10. 
0.0091 0.0097 0.0102 0.0157 0.0213 0. 
2.1166 2.0586 2.0014 1.4696 1.0123 Q. 
3.9520 4.0403 4.1245 4.7338 4.9209 4. 
-4.0001 -4.1715 -4.3355 -5.5627 -6.0402 =5. 
-0.5203 -0.4271 -0.3385 0.2995 0.4859 0. 
0 0 0 0 0 
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oh hOoOwonooco 





The following is a partial listing of the output of a 


MATLAB function that gives the percent error in the desired 
element value and the value calculated from the coefficients 


at specific geometries. 


first five places represent, 
places represent the element values. 


The data is arranged in columns, the 


a,b,W,d,T and the last five 
As discussed above only 


the inductor shows any significant deviation from the desired 
values. 


element_error = 


Columns 1 


900. 
.0000 
400. 


400 


10. 
. 3010 
-0. 

0. 
-0. 
.0000 
Q. 


~17 


0 


Columns 8 


900. 
400. 
400. 


0000 


0000 
0 
0000 


0000 
0000 
0000 


0000 


0000 
0000 
0000 

0 


.0000 
- 7242 
.0000 
-9000 
-0000 
-0000 
-0000 


through 7 


900.0000 900. 
400.0000 400. 


400.0000 400. 
0 
20.0000 30. 
-13.6287 -10. 
-0.0000 -0. 

0.0000 QO. 
-0.0000 ~0. 

0.0000 ~0. 

0.0000 0. 
through 14 


900.0000 900. 
400.0000 400. 


400.0000 400. 
0 
90.0000 100. 
~7.8888 -8. 
-0.0000 aes 
~0.0000 -0. 
~0.0000 =O. 
0.0000 QO. 
-0.0000 =O: 
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Columns 15 through 21 


900.0000 
400.0000 
400.0000 
10.0000 
10.0000 
-5.0626 
0 

0.0000 
-0.0000 
-0.0000 
0.0000 


. 0000 
.0000 
.0000 
.0000 
.0000 
.0575 
.0000 
0000 
.0000 
0000 
.0000 


Columns 22 through 28 


900.0000 
400.0000 
400.0000 
10.0000 
80.0000 
12.2612 
-0.0000 
-0.0000 
-0.0000 
0.0000 
-0.0000 


.0000 
.0000 
.0000 
.0000 
.0000 
.9748 
.0000 
.0000 
0000 
.0000 
.0000 


Columns 29 through 


900.0000 
400.0000 
400.0000 
30.0000 
10.0000 
-4.9390 


-0000 
.0000 
.0000 
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.0000 
.0000 
.0000 
.0000 
.0000 
. 3605 
.0000 
.0000 
.0000 
.0000 
.0000 


.0000 
.0000 
.0000 
.0000 
.0000 
. 8648 
.0000 
.0000 
.0000 
.0000 
.0000 


.0000 
.0000 
.0000 
.0000 
.0000 
.3735 
.0000 
.0000 
.0000 
.0000 
.0000 


.0000 
.0000 
.0000 
.0000 
.0000 
.0122 
.0000 
.0000 
.0000 
.0000 
.0000 


.0000 
.0000 
0000 
.0000 
.0000 
.5673 
.0000 
.0000 
.0000 
-0000 
.0000 


.0000 
.0000 
.0000 
.0000 
.00C0 
. 7318 
.0000 
.0000 
.0000 
.0000 
.0000 


.0000 
.0000 
.0000 
.0000 
.0000 
. 2380 
.0000 
.0000 
.0000 
.0000 
.0000 


.0000 
.0000 
.0000 
.0000 
.0000 
.4440 
.0000 
.0000 
.0000 
.0000 
0000 


.0000 
.0000 
.0000 
-0000 
.0000 
.7550 
.0000 
.0000 
.0000 
.0000 
.0000 





Columns 36 through 42 


.0000 
-9000 
.0000 
.0000 
.0000 
.7284 
.0000 
.0000 


.0000 
.0000 
.0000 


.0000 
.0000 
.0000 
.0000 
.0000 
- 4490 
.0000 
.0000 





.0000 
.0000 
.0000 
.0000 
.0000 
. 1606 
.0000 
- 0000 


-0000 
.0000 
-0000 


900.0000 
400.0000 
400.0000 
30.0000 
200.0000 
—-3.8652 
-0.0000 
-0.0000 


0.0000 


0.0000 
-0.0000 
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.0000 
.0000 
.0000 
.0000 
.0000 
-0880 
.0000 


0000 


.0000 
.0000 
.0000 





Columns 1 


900.0000 
400.0000 
400.0000 
0 
10.0000 
14.5784 
0.0047 
1; 5939 
2.5056 
1.5626 
-3.2824 
30.7504 
4.9488 


The following is a listing of the output of a function 
that iterates through the STRIP data and gives the resulting 
Maximum percent error for that geometry with the model using 
element values generated using coefficients. 

The columns are arranged in the following way: 


a 
b 
Ww 
d 
T 
L 


inductor in nH 


C capacitor in pF 


0 order coefficient for the turns ratio 


lst order coefficient for the turns ratio 
2nd order coefficient for the turns ratio 
3rd order coefficient for the turns ratio 
% error phase 

% error magnitude 


through 7 


900.0000 900. 
400.0000 400. 
400.0000 400. 
0 
40.0000 70. 


14.9579 15. 
0.0054 0 
1.2798 1 
3.7159 4 

-1.2922 =33 

-1.5344 Q. 

23.6768 23. 
6.5312 9. 

through 14 


10.0000 10 
200.0000 500 
12.3204 ~8 
0.0079 0 
0.6374 0 


0000 
0000 
0000 

0 


0000 
0325 


.0063 
0331 
.8164 


8981 
0636 
2986 
8957 


.0000 
-0000 
.0000 
-0000 
.0000 
- 4988 
.0196 
- 8890 
- 9289 
-7254 
-4575 
.1516 
.9739 


900. 
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0000 


.0000 
.0000 
-0000 
.0000 
.7420 
.0048 
.4271 
-4819 
. 1668 
-4472 
.6509 
. 7641 


.0000 
.0000 
-0000 
.0000 
.0000 
.2218 
.0045 
.3245 
-7701 
. 3334 
.6914 
.8377 
- 7996 


.0000 
.0000 
.0000 
.0000 
.0000 
- 8628 
-0051 
-2465 
.0648 
7039 
. 3651 
. 6462 
.8559 


.0000 
.0000 
. 0000 
-0000 
. 0000 
-5792 
.0081 
.1507 
.1144 
.0978 
.1394 
.2073 
. 7834 





Columns 15 through 21 


185 


-0000 
.0000 
.0000 
.0000 
.0000 
- 6882 
.0050 
. 6733 
. 6350 
.0731 
. 9066 
.0244 
61. 


7791 


111 


.0000 
.0000 
.0000 
.0000 
.-0000 
. 3058 
.0048 
. 6406 
-0694 
-1116 
. 6769 
7432 
65. 


7842 


Columns 22 through 


.0000 
.0000 
.0000 
.0000 
.0000 
.2627 
.0063 
.8584 
-4892 
. 2833 
-0488 
.7916 
-1136 


.0000 
.0000 
.0000 
.0000 
.0000 
. 3093 
.0094 
-5995 


263:74 


hrs 


.1580 
-2808 
- 4038 
-1170 


Columns 29 through 35 


.0000 
.0000 
.0000 
.0000 
.0000 
-4258 
.0051 
.6570 
-9855 
. 5888 
.7114 
.9502 
. 7010 


.0000 
-0000 
-0000 
-0000 
.0000 
-4140 
.0051 
. 6234 
-0660 
-7776 
- 5832 
.6716 
-6612 


.0000 
.0000 
.0000 
.0000 
.0000 
-5639 
.0048 
.5819 
.1867 


6455 


. 3325 
7305 
- 6520 


.0000 
.0000 
.0000 
.0000 
.0000 
.8442 
.0064 
.4179 
- 4468 
.0304 
.0538 
.0123 
1120 


.0000 
.0000 
-0000 
-0000 
-0000 
.8549 
-0054 
.1991 
.1763 
.3961 
-2304 
-6615 
. 4526 
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.0000 
-0000 
.0000 
-0000 
-0000 
3139 
0136 
.4235 
.7394 
-4551 
-9145 
-0340 
.6357 


0000 


.0000 
.0000 
.0000 
.0000 
.1552 
.0080 
.2349 


7625 


.6347 
.7201 
.2764 
-5835 


0000 


.0000 


0000 


.0000 
.0000 
sod 
.0056 
.0682 
. 6334 
-4885 
.0265 
.0107 
.9895 


.0000 
-0000 
.0000 
.0000 
.0006 
.7979 
.0057 
~9275 
. 7895 
.8210 
1922 
1165 
. 3909 


.0000 
0000 
.0000 
.0000 
.0000 
.3745 
0213 
0123 
. 9209 
.0402 
- 4859 
.7164 
1735 


.0000 
.0000 
-0000 
.0000 
-0000 
.2001 
.0043 
.3707 
9955 
- 8458 
.8102 
.2965 
.9708 
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